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It is common practice to estimate the errors on the angular power spectrum which could be obtained 
by an experiment with a given angular resolution and noise level. Several authors have also addressed 
the question of foreground subtraction using multi-frequency observations. In such observations the 
angular resolution of the different frequency channels is rarely the same. In this report we point out 
how the "effective" beam size and noise level change with £ in this case, and give an expression for 
the error on the angular power spectrum as a function of I. 

PACS numbers: 98.70.Vc, 98.80.Cq 
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I. INTRODUCTION 

There are many experiments planned which hope 
to reap the rich harvest of information available in 
the anisotropies of the Cosmic Microwave Background 
(CMB) |]. For the purposes of estimating cosmologi- 
cal parameters or constraining models it is important to 
know how well an experiment can constrain the angular 
power spectrum of CMB anisotropies. 

It is common practice to estimate the errors on the an- 
gular power spectrum which could be obtained by an ex- 
periment with a given angular resolution and noise level. 
Such an estimate can give us insight into what regions 
of the angular power spectrum one could constrain, and 
what are the limiting factors in an experiment designed 
e.g. to constrain cosmological parameters. In this report 
we point out some simple limiting cases which allow one 
to gain intuition about the effect of multiple observing 
frequencies with differing beam sizes in the presence of 
(somewhat idealized) foregrounds. 

If we write the angular power spectrum of the 
anisotropy as Cg and the noise power spectrum as Ng 
(typically Ng — 47r/ s i cy cr 2 /A / pi x , where / s k y is the fraction 
of sky covered, a is the rms pixel noise and N p i x is the 
number of pixels) then 



SCg 



c +^ 
(2£+l)/ sky V * Wg 



(1) 



where Wg is the window function of the experiment and 
we have assumed that the noise is gaussian. For a gaus- 
sian beam of width 9b, Wg = exp(— £ 2 9l). Eq. (|J) as- 
sumes that the only sources of variance in an experiment 
are the anisotropy signal and the receiver noise, i.e. it ne- 
glects foregrounds. We have also implicitly assumed that 
SCg is a gaussian error, thus we imagine we are work- 
ing at reasonable high-£ and will bin our power spectrum 
estimates into finite width bins in I. Eq. (fy) has been 



widely used to estimate how well upcoming satellite ex- 
periments could constrain cosmological parameters |1 . 

Several authors have also addressed the question of 
foreground subtraction using multi-frequency observa- 
tions ||-|(| • With multi- frequency observations it appears 
possible to separate the desired anisotropy signal from 
the foregrounds with encouraging precision. The authors 
of |^ have shown that one can regard foreground sub- 
traction as an enhancement of the noise in a foreground 
free experiment. The noise enhancement factor has been 
called the foreground degradation factor (FDF) by Dodel- 
son jj] who gave a simple expression for it. A method of 
foreground subtraction using both frequency and spatial 
information was proposed in Ref. ^j. Various methods 
of foreground subtraction have been compared recently 
by Tegmark @. 

Note that the "noise term" in Eq. (Q) depends expo- 
nentially on I once t > O^ 1 . In multi-frequency obser- 
vations the angular resolution of the different frequency 
channels is rarely the same, leading to the question of 
which beam size to use. In this report we discuss a sim- 
ple heuristic "effective" beam size and noise level. Since 
how well we can take out foregrounds depends on the an- 
gular structure of the foregrounds, the "noise" will be a 
function of the foreground angular power spectrum. 

The formalism is directly comparable to that in [^|J6| in 
that we use a minimum variance estimator of the CMB 
anisotropy power spectrum. An extension to "real world 
foregrounds" p] is straightforward, but the purpose of 
this report is to gain intuition through simple examples 
so we do not pursue this line of development. 



II. FORMALISM 

As noted by Tegmark || , one can consider foregrounds 
as an additional noise component which is highly corre- 
lated among different frequency channels. If we label the 
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frequency channels by a greek subscript we can define a 
noise correlation matrix 
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which has contributions from the pixel noise (assumed 
uncorrelated between channels here for simplicity) and 
foregrounds labelled by superscript i. Here (f} a fla) is 
the correlation matrix of foreground i in channels a and 
/3, with the 4n/Npi x inserted for later convenience. If the 
foregrounds are 100% correlated between the channels 
then (fiJip) = fLftf} where fj a is the rms intensity as 
a function of frequency. We shall assume this case from 
now on, but see Ref. M. 

It is straightforward to derive the minimum variance 
estimator of Ce, as a linear combination of measure- 
ments at different frequencies. If we have measured mul- 
tipole moments af m at frequency a then we write the 
estimate of the CMB component as 0e m = J2 a F a a" m . 
Imagine that we can write the observed signal af m = 

1/2 

timWfa + nemo, where tg m is the cosmological signal of 
interest (the same in each frequency channel) and ng ma is 
the sum of the noise and foregrounds (i.e. the non-CMB 
components). We now minimize the variance of our es- 
timator 9g m minus the "real" underlying sky (te m ) with 
respect to the weighting matrix F a . We find that our 
minimum variance Ce estimator is (averaging over m) 
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(3) 
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where Of is the average of 6f m over 
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and we have left the ^-dependence of F a implicit for no- 
tational convenience. The vector F a projects out the ^th 
CMB multipole moment from the signal in each channel 
in a minimum variance sense. 

If we assume that this estimator is Gaussian then we 
can replace Eq. (|l|) by 
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where we have set / s k y = 1 for simplicity (the scaling 
with / s ky is given in Eq. (jl])). Note that in the limit of 
one frequency channel and no foregrounds, the sums over 
a and (3 are trivial, the noise term is independent of F 
and we recover Eq. (Q). 

Eq. (^]) is the general result, and we consider several 
examples to gain intuition in the next section. 



III. EXAMPLES 

Let us consider various limits of Eq. (§). We will focus 
on the noise term, since the part of the error proportional 
to Ce simply reflects cosmic plus sample variance . In 
the signal dominated limit [Ct ^> Ng) F a is just the 
inverse square root of W( a and the noise term in Eq. (|^) 
reduces to 
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where We = J2 a Wa- Thus, in the absence of fore- 
grounds the noise is the sum of the noises in those chan- 
nels able to resolve features with multipole number £, 
divided by the number of such channels squared. The 
addition of foregrounds increases the variance, but by 
assumption both contributions to SCi are sub-dominant. 
We expect this regime to occur at low-£, where we are 
cosmic and sample variance limited. 

Now let us consider the opposite limit. First imagine 
there are no foregrounds. In this limit (Ne ^> Cg) the 
noise term in Eq. (|^) reduces to 
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which would be the obvious way to combine the chan- 
nels to obtain the effective noise: recall the error on a 
weighted mean is a~ 2 — <r^ 2 . 

Now we can enhance this last example by the addition 
of foregrounds. For simplicity imagine a two channel ex- 
periment with 1 foreground fe a in addition to uncorre- 
lated noise a a , with a = 1,2. For dcfinitcncss imagine 
that channel 1 is a low frequency channel with a "large" 
beam, while channel 2 is a high frequency channel with a 
"small" beam. Some trivial matrix algebra allows us to 
write [Ne/We) c g as 
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iV pix u\W n + a\Wex + [fn - fa] 2 WnW a 



(8) 



which reduces to Eq. (0) as f a — > 0. 

To bring out the essential details let us take = er 2 
and ignore the different resolutions by working at low 
enough £ that We a — ► 1, 
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The term in parentheses is the increase in the noise 
over the one channel result. As we decrease a 2 , this goes 
from | (the case of no foregrounds: co-add the chan- 
nels) to (1 + x 2 )(l — x)~ 2 where x = jajjti. For this 
case, the increase in the noise is just the FDF defined 
by Dodelson p}, written in a slightly different notation. 
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The point x = 1 is when the frequency dependence of 
the foreground is the same as the CMB. Generically the 
minimum variance estimator Eq. (|J) does better than the 
FDF would indicate, as has been pointed out in Ref. ||. 

Now let us reinstate the window functions, but imagine 
we are at intermediate £, where Wn <C 1. In this regime 
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So the appropriate resolution is that of channel 2, and 
the noise is the channel 2 noise plus a contribution from 
the foreground at the higher frequency. Note that the 
properties of channel 1 do not enter the expression, as 
expected. The foreground contribution declines as one 
approaches the resolution of channel 2, so that asymp- 
totically the error is just 47rcr|/iVpi x W / £2 • 

Finally note that if our foregrounds have a "steep" 
power spectrum, i.e. fall rapidly with £, then at high-£ we 
have the limit fg a — > and reproduce Eq. (^|). At low-£ 
the foreground should dominate over the noise which is 
the case discussed below Eq. (|^). This shows that for 
foregrounds without a lot of small-scale structure there 
is little effect at high-i? from removing the foregrounds (in 
this simplified example where the foreground properties 
are assumed known exactly). 
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IV. SUMMARY 



Let us summarize our results with I increasing from 
the signal- to the noise-dominated limits. Under our 
assumptions, the error on the angular power spectrum 
starts dominated by cosmic and sample variance at low- 
£. Moving to higher £, if the noise or foregrounds start to 
dominate on scales which all channels can resolve, the ef- 
fect of a foreground looks like an increase in the noise by 
a factor which is always less than Dodelson's FDF j| . If 
the noise or foregrounds dominate on scales smaller than 
the size of the widest beam, but larger than the size of 
the smallest beam, the noise is that of the channel with 
the smaller beam size, increased by the variance of the 
foreground at that frequency. At £ > B^ 1 for even the 
highest frequency channel the error is simply that ob- 
tained from the noise in the highest frequency channel, 
regardless of foregrounds. If the foregrounds are always 
negligible compared to the noise, the appropriate noise 
level is the weighted sum of the noises in each channel 
which can resolve the given £, as shown in Eq. (R). 
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